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Abstract 
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m ; 

£— : ■ We re-analyse the long distance contributions to the process b — ► sl + l~. Full (j 2 - 

qv . behavior of the vector meson dominance amplitude is used together with the effect 

of Terasaki suppression, and comparisons with the previous results are given. We 

show that the interference between short- and long distance contributions makes it 

Q : 

_C ■ difficult to extract the short distance information from the dominant long distance 

> ■ 

background, either in the dilepton invariant mass distribution or in the single lepton 

><■ 

energy spectrum. 
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Rare decays through the flavor changing b — > s transitions provide good test of the 
standard model (SM), and are expected to give signals of new physics [[I]]. The branching 
ratio of the process b — > sj, which has already been measured by the CLEO collaboration 0, 
is within the SM predictions. Unlike the decay b — ► 57, the process b — > sl + l~(l = e or 
/i) is expected to be dominated by long distance contributions through the mechanism of 
vector meson dominance(VMD)P]. However, it was usually believed that the long distance 
(resonance) contributions arise only in some particular region of the invariant mass spectrum 
of the dilepton pair||], since the involved resonance ip(ip t ) peak is very sharp. Detailed 
calculation || shows that there exists significant interference between the short and long 
distance contributions, which leaves only a small portion of kinematic region at low dilepton 
invariant mass where the interference effect by the resonances is small. The energy spectrum 
of single lepton has also been given in || where a window of nearly pure short distance 
information is found. 

In the present work, we will re-examine both the dilepton invariant mass and the single 
lepton energy spectrums using an alternative treatment of the VMD amplitude. In the pre- 
vious analysis of the cascade decays b — » sip(tp') and ip{ip') —> l + l~ , an effective description 
is made for the later electromagnetic transition ip(ip t ) — > / + /~[[|, f|. The dependence of the 
VMD amplitude on the square of the dilepton invariant mass, q 2 , is approximated by that 
of the resonance mass vnL or mi in the denominator of the photon propagator. This ap- 
proximation is only valid near the resonance region, and consequently, the previous analysis 
are not complete in the whole phase spaceQ 

Let us start with the short distance contributions to b —>■ sl + l~ with I = e or jjl. The 
short distance contributions come from box, Z and photon penguin diagrams. The QCD 
corrected effective Hamiltonian in SM is HI: 
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lr This has also been noted in |6| in the discussion of exclusive decays of B meson, and in J/j in the charm 
quark sector. 



with P L = (1 — 7 5 )/2, Pr = (l + 7 5 )/2, and q = P1++P1- is the invariant mass of the dilepton. 
The analytic Wilson coefficients Cj (fi), C$ (//), and Cio(/x) are given in Ref.||. Under 
the leading logarithmic approximation, we get the numerical results at \x = m b as: 

(% ff = -0.315, C 10 = -4.642, (2) 

and to the next-to-leading order, 

C 9 e// = 4.227 + 0.124 u(s) + 0.359 g(m c /m bl s)+ 0.034 g{l,s) + 0.033 #(0,s), (3) 

where s = g 2 /m 2 . The function w(s) and g(z,s) can be found in ref.fpfl. Here for numerical 
evaluation, we use m top = 176GeV§, m b = 4.8GeV, m c = 1.4GeV, A QCD = 225MeV 0. 



By normalizing to the semileptonic rate, the strong dependence on the b-quark mass 
cancels out. The differential decay rate dT{B -^ XJ + l~)/ds, where s = (pi+ + Pi-) 2 /m%, is 
given by 



d T(B^X s l + l-) = —^£ 7 ^(l-s) 2 \(l + 2s)(\C e 9 ff \ 2 + C 2 w 
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where f(m c /mb) is the phase space factor: 

f(x) = 1 - 8s 2 + 8x 6 -x 8 - 24x 4 In x. 



If we take the experimental result Br(B — > X c eu) = 10.8%[10|, the differential decay rate of 



B — > X s ji + n is found, which is depicted in Fig.l as the dash-dotted line. 

In addition, there are also long distance resonance contributions from cc state. There 



are six known resonances in the cc system that can contribute to this decay mode |TT 
The lowest two, ip and ip' , were considered in the previous analyses [f|, ||. Here we also 
consider the same two resonances. The higher resonances will also contribute, but they are 
less important in our case of discussing the uncertainties, as will be shown later. 



Applying the VMD mechanism, the long distance contribution is through b —>■ sip, and 
ip — > 7 — > l + l~ , where the resonance can also be ip 1 . These give the effective Lagrangian 

Cres = ^ftf) ^V*K(*fPLb)hj + W - VO, (5) 

3q 2 (q 2 -m^ + im^T^) y/2ir 

where 02 = C\ + C2/3 is a QCD corrected coefficient of the four quark operators. Below we 
will use the phenomenological value 02 = 0.24 which comes from fitting the data of B meson 
decays []TJ]. Note that the expression (|) differs from the previous ones|3|, Q by keeping the 
photon propagator as —ig^ u /q 2 instead of —ig^/wL or —ig^/rrAi. Thus it holds in the 
whole kinematic region. 

The effective coupling of a vector meson gv(q 2 ){V — ip , ip' ) is defined by 

<0\c lfl c\V(q)>=g v (q 2 )el, (6) 

where eY is the polarization vector of the vector meson V. On the mass-shell of the vector 
meson, gv{q 2 ) is replaced by the decay constant 5v( m y) which can be obtained from the 
leptonic width of the vector meson: 

T{v^tn = ^r/vi<). (7) 

The structure of eqn.(^) is the same as that of the operator 9 . It is convenient to include 
the resonance contribution in eqn.(3) by simply making the replacement 

ci" - <%" = cr + p!^m + (v _ n (8) 

3q 2 (q 2 -m^ + im^T^) 

Assuming a constant coupling g"^(q 2 ) = g%{mV) as done in [0, 0, the numerical result is 
given in Fig.l as the dashed line. It is easy to see that, this spectrum is enhanced in the low 
q 2 region due to the explicit inclusion of the photon propagator. From Fig.l, we can also 
expect that higher resonances other than ip or ip' contribute mainly in the region 0.6 < s < 1 
where we are not interested, since near this tail of the spectrum no useful short distance 
information is expected to emerge. 



The assumption made above on the constant coupling of g^(q 2 ) can be improved by ac- 
counting for the mechanism of Terasaki suppression for the ip — 7 conversion ||13|| . In the 
framework of VMD, the data on photoproduction of if) indicates a large suppression of g^(0) 
compare to p^(m3,)|T^]. This has been confirmed in |DJ by constraining the dominant long 
distance contribution to s — > dj using the present upper bound on the Q~ — > 2~7 de- 
cay rate. As a result, it can be concluded that this suppression results in a much smaller 



long distance contribution to b — > 37 transition |14J. Now we use a momentum dependent 
9v(q 2 )(V — ip j i>') i n ^res, which is used in [T^ to obtain a reduced resonance to non- 
resonance interference where a broader region of invariant mass spectrum sensitive to short 
distance physics is claimed. 

The momentum dependence of gy(q 2 )(V = ip , ip') derived using a dispersion relation JIB | 
is 

gv(q 2 )=gv(0)(l + ^[dv-h(q 2 )]Y (9) 

where Cy, = 0.54 , c^t = 0.77 and d$ = dy = 0.043. h(q 2 ) is defined by 



h ^ -v^-r\- i - 2 f + ^ + 2r V; - ' tan " 1 jh \ (10) 



with r = q 2 /my for < q 2 < rriy. As a result, eqn.@, which is valid for < q 2 < my, is an 
interpolation of gy from the photoproduction experimental data on gy(0) to gy(m v ) from 
the leptonic width based on quark- loop diagram. We assume gy(q 2 ) = gy{m v ) for q 2 > m v 
mainly due to the fact that the behavior of the ip — 7 conversion strength is not clear in this 
region, and is not important in our case (see below). 

Applying Terasaki's formula (0) for the q 2 dependence of gy(q 2 ), the differential decay 
rate of b —>■ sl + l~ receives suppression in low q 2 region. However, there is still significant 
interference between the resonance and the short distance contributions, due to the factor 
1/q 2 coming from the propagator of the virtual photon. This is also shown in Fig.l. 

Now we turn to the energy spectrum of single lepton. The integration over q 2 is compli- 
cated, since many functions here involve s. We simply give the numerical results in Fig.2 for 



I = /z(see also |J). One can see that, if the 1/g 2 behavior is replaced by 1/rri^ (or l/mL) 
everywhere, there is almost no contribution from the if), if)' resonance when f3 < 0.2. This 
result is what has been arrived in Ref.||. If, however, this 1/q 2 is retained, there are also 
contributions from the resonances even in the low (3 region and consequently, the resonance 
background is still serious. Including the effect of Terasaki suppression of gy(q 2 ) in the low 
q 2 region, the result is also shown in Fig. 2 where the resonance background is only half 
reduced. 

From both Fig.l and Fig. 2, we can observe that the long distance VMD contributions 
to the process b — > sl + l~ are large, if an alternative treatment of the electromagnetic sub- 
process is performed. It can also be seen that the single lepton energy spectrum is almost 
useless in extraction of any short distance information from the resonance background. The 
total branching ratio of b — ► sl + l~ turn out to be 3.6 x 10~ 4 or 4.9 x 10 -4 with the effect of 
Terasaki suppression included or not. 

We have treated the resonance contribution from if), if}' to b — > sl + l~ alternatively without 
using the effective description of the electromagnetic sub-process if)(if)') — > 7/ + / _ . The long 
distance contributions are found to be significant, especially in the single lepton energy 
spectrum. Concerning other higher resonances, and also other uncertainties existed in this 
decay mode||, we conclude that it is difficult to extract short distance information, which 
is sensitive to new physics, from the dominant long distance contributions. 
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Figure Captions 

Fig.l The differential decay rate via s = m 2 /ml. The dash-dotted and the dotted lines 
correspond to the results without resonance contribution and with resonance contribution 
included as [§, ||, respectively. The dashed and solid lines are the results with resonance 
included with the treatment of the electromagnetic sub-processes using eqn. (5), while 
constant coupling gv(q 2 ) = fiv( m v) an d © ar e used, respectively. 

Fig. 2 Same as Fig.l for the differential decay rate via (3 = E~ jm^. 
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